1. IntroductionTransmutation by neutrons, produced by high flux accelerated protons bombarding a heavy metal target, is the primary method to deal with long-lived fission products (LLFPs) with large neutron capture cross-sections. However, for nuclides, whose neutron capture cross sections are small, this method is not very efficient. In addition, (n, γ) reactions may transmute stable nuclides into long-lived radiotoxic nuclides. For example, 133Cs is stable, of which the (n, γ) transmutation product will be 134Cs followed by 135Cs. They are unstable and radiotoxic with long half-lives of 2.06 years and 2.3 million years, respectively.
Photo-transmutation can be a good supplement to neutron transmutation.[1] There are several methods to obtain γ rays: positron annihilation in flight, bremsstrahlung, laser-Compton scattering (LCS), and nuclear resonance reactions. γ-ray beams obtained by positron annihilation in flight is of low intensity, thus it is not appropriate for transmutations. Transmutations of 129I based on a laser-induced bremsstrahlung (LIB) source have been already performed experimentally in Germany,[2–4] UK,[5] and Japan.[6] Theoretical investigations have been done on LLFPs such as 126Sn,[7,8] 107Pd,[9] 93Zr,[10] 90Sr,[11] 135Cs,[12–14] and 137Ca.[15] An LCS γ-ray source can produce γ-rays in a wide energy range, which can be utilized to transmute LLFPs. The Upgraded HIγS, operating since 2009, can produce up to 3×109 photons/s. The energy of photons can be adjusted in 1 MeV–100 MeV. The energy of electron beam is 0.24 GeV–1.2 GeV and the laser wavelength is
. The rate of transmutation reaction of 135Cs by the Upgraded HIγS is evaluated in Ref. [16] to reach 1011 per hour. Studies[17,18] show that the LCS-based photo-transmutation is more efficient than LIB. Many experimental[19,20] and theoretical[14,16–18,21] studies have been done using LCS γ-ray source to transmute LLFPs.
In the present paper, γ rays emitted from resonance reactions induced by proton is utilized to transmute LLFPs. Nuclear resonance reaction can produce mono-energetic γ rays, which can be used to transmute LLFPs. To do so, the energy of incident protons is supposed to be hundreds of keV to several MeV. Figure 1 shows high power proton beams worldwide (non exhaustive) crossing this region.[22–30] The mean beam current of high power proton beams that have already been constructed can reach 100 mA.
2. Proton-induced resonance γ-ray source for transmutationThe excited levels of the compound nucleus are discrete and widely spaced when the bombarding proton energies are low. The so-called Breit–Wigner formula, as shown below, gives a good description for the cross section of a single resonance level,
where
σ(
Ep) is the (p,
γ) reaction cross section of the protons whose energy is
Ep in laboratory frame,
Ecm,p is the corresponding energy of protons in center-of-mass frame,
,
E0 is the resonance energy.
ma,
mA are the mass number of the incident proton and target nucleus.
is the statistical factor,
IC,
Ia,
IA are angular momentum of the compound nucleus, the proton and target nucleus, respectively.
λp is the de Broglie wavelength of relative motion, Γ
T is the full width at half maximum (FWHM) of the compound nucleus level. Γ
γ is the partial width for the gamma ray transition to the ground state and Γ
p the partial width for proton emission. Thick target
γ photon yields of some nuclides are shown in Table
1.
[31–33]Table 1.
Table 1.
| Table 1.
γ yields of some narrow resonance reactions.
. |
The GDR peak energy of most LLFPs is about 15 MeV, therefore it would be better to use the reaction 7Li
8Be to obtain higher transmutation rates. According to the data listed in Table 1, the yield of 14.8 MeV (33%) and 17.6 MeV (67%) γ rays emitted from 7Li
when the energy of incident protons is 441.4 keV is 1.90×10−8 per proton in total.[31,32] The yield of 9.17-MeV γ-ray emitted from 13C
14N is 4.7×10−9.[33]
Different compound nuclides have different energy levels ranging from several hundreds of keV to tens of MeV, which cover most Giant Dipole Resonance (GDR) regions. The yields of γ-radiation produced by proton beams in narrow resonance region follow the equation:
where
is the number of resonance energy
γ rays produced per second in the solid angle
,
is the intensity distribution against the incident proton energy, which is a function of the depth of the target l,
is the differential cross-section when the proton energy is between
E and
and the emitted
γ rays is in the solid angle
to
,
n1 is the number density of the target atom in the converter.
l1 is the thickness of the converter,
E1 and
E2 are the upper and lower bounds of the incident proton energy.
Next, using the γ-radiation produced above, we obtain the transmutation rate without considering secondary process:
where
nt is the transmutation rate,
is the number of
γ photons produced per second in the solid angle corresponding to the transmutation target,
is the attenuation coefficient of the target which is relevant to the
γ photon energy
E0,
l2 is the thickness of the target,
is the cross-section of the photo-nuclear reactions when the energy of incident
γ is
E0,
n2 is the number of nuclei to be transmuted. However, this is a complex process.
γ photons can produce electrons and positrons through pair production process, these electrons and positrons produce photons again by bremsstrahlung, which can also add to transmutation rate. It is better to track all the particles with a Monte Carlo simulation.
From Eqs. (2) and (3), we can increase the transmutation rate nt in the following ways: 1) higher proton beam intensity Ip; 2) reactions with higher cross-section of proton induced γ-rays; 3) larger solid angle of the transmutation target.
3.
197Au
196Au experimentFirst, to ensure the accuracy of the calculation, the experiment of 197Au
196Au was carried out in the 2×1.7-MV tandem accelerator of China Institute of Atomic Energy (CIAE).[33] Two gold targets of the same size are irradiated. Figure 2 shows the transmutation process. The γ photons are produced by 13C
14N reaction when the resonant energy of incident protons is 1.747 MeV. 13C is evaporated onto a ϕ10×2 mm gold target up to a mass thickness of
. Bombarded by the 9.17-MeV γ photons, 197Au can be transmuted into 196Au through
reaction. However, 196Au is unstable, it decays to 196Pt through electron capture process and to 196Hg through β decay.
The first sample is radiated for 6 hours and the second for 5.5 hours. Then, the radioactivity of the irradiated gold samples is measured in a low-background lead chamber with anti-Compton method used to reduce the background radioactivity. The detected counts obey the exponential law given by
where
P is the generation rate during irradiation and
λ the decay constant of of
196Au.
is the irradiation time,
is the cooling time, and
is the measuring time. By measuring the decay counts
Nc, the generation rate
P can be derived from Eq. (
4), thus the transmutation counts
N0 can be obtained from
. The spectrum is presented in Fig.
3. Peaks at 355.73 keV (I
γ=87%) and 333.03 keV (I
γ=22.9%) originates from
196Au. The result is shown in Table
2. In Table
2, the second column is the irradiation time, the third column shows the counts of peak at 355.73 keV, and the numbers of (n,
γ) in
197Au, calculated from experimental data, with self-absorption of the gold targets taken into consideration, are listed in the fourth column.
Table 2.
Table 2.
Table 2.
Transmutation numbers of different samples.
.
Sample |
Tirradiation/s |
NDet
|
Ntr/106
|
1 |
21600 |
1566 ± 49 |
1.01 ± 0.03 |
2 |
19800 |
6610 ± 100 |
1.09 ± 0.02 |
Average |
20700 |
– |
1.05 ± 0.03 |
| Table 2.
Transmutation numbers of different samples.
. |
The thick target γ photon yield of 13C
was measured with an HPGe detector (GXM35P4-70). The maximum yield of 9.17-MeV γ photons is (4.7±0.4)×10−9/proton.[33] As the average current of the proton beam during the 197Au transmutation experiment is
, the number of incident protons during 20700-s irradiation time is 1.04×1018, producing (4.9±0.4)×109 γ photons.
In the Geant4 simulation, 4.9×109 γ photons are assumed to be emitted in 4π direction from a point source. The point source is on the center of the front surface of a ϕ10×2 mm gold target, with zero distance. The number of transmutation is (6.6±0.5)×105, which is 37% less than the experimental results but with no deviation in magnitude.
4. Geant4 simulation of photo-transmutation yieldsIn this section, we simulate photo-transmutation yields of nuclear waste like 129I. The photonuclear cross sections of 129I extracted from Geant4 packages are shown in Fig. 4.
Photoelectric effect, Compton effect, pair production, and photo-nuclear reactions are the principal mechanisms by which γ rays interact with matter. Figure 5 shows cross section of each mechanism when γ rays interact with 129I.
Assuming the intensity of the proton accelerator is 10 mA at 441 keV, according to the yields of 441-keV narrow resonance of reaction 7Li
8Be in Table 1, γ-ray yield of 14.8 MeV and 17.6 MeV can reach 1.2×109 per second, and the γ-ray source, produced by the reaction 7Li
8Be, is assumed to be point shape and isotropic in 4π direction, of which 33% are 14.8 MeV and 67% are 17.6 MeV. Next, the 129I target, with width and height are 2 cm, is placed 5 cm away from the point source. The 129I
128I yields against target depth are shown in Fig. 6, the maximum yield can reach 5.7×10−5 per γ. Taking the γ-ray yield into consideration, there are 6.8×104 129I
128I reactions per second.
Regarding transmutation of other LLFPs based on 7Li
8Be 441-keV resonance γ-radiation in the same condition as mentioned earlier but the point source with zero distance is on the front surface center of the targets and the targets are as thick as 10 cm. The rates of transmutation reactions are calculated and shown in Table 3.
Table 3.
Table 3.
Table 3.
Transmutation rates of different LLFPs.
.
Reactions/hour |
129I |
135Cs |
126Sn |
93Zr |
107Pd |
137Ca |
90Sr |
|
4.7×109
|
2.1×109
|
2.8×109
|
8.6×1010
|
9.7×1010
|
1.8×109
|
2.8×109
|
|
4.7×109
|
2.2×109
|
7.9×109
|
7.9×1010
|
5.4×1010
|
2.5×109
|
6.1×109
|
Total |
9.4×109
|
4.3×109
|
1.1×1010
|
1.7×1011
|
1.5×1011
|
4.3×109
|
8.9×109
|
| Table 3.
Transmutation rates of different LLFPs.
. |
The transmutation rate of 135Cs is 109 per hour is shown in this chart in total. According to Ref. [16], among LCS facilities constructed around the world, the best transmutation ability is of the Upgraded HIγS, which can reach 1011 per hour, the HIγS is 109 per hour, and the NewSUBARU is 108 per hour.
5. ConclusionConsidering a 10-mA proton accelerator, nuclear waste transmutation driven by 7Li
8Be 441-keV narrow resonance γ-rays have been investigated through Geant4 simulation. The maximum transmutation yield of 129I is 9.4×109 per hour. The transmutation yield of 135Cs is obtained to be the same order as HIγS of 109, which is less than that of the upgraded HIγS by one order of magnitude but greater than that of the NewSUBARU by two orders of magnitude. In the future, we will improve the transmutation rates by improving the current of the proton accelerator, searching for higher yields of γ-radiation in GDR region, designing better target structure, and so on. Moreover, when compared with the good direction of LCS γ source, the γ-radiation emitted from narrow resonance is in
directions, which is appropriate for large areas of LLFP radiation.
Acknowledgment
The authors would like to thank Prof. Bing Guo, Prof. Hong-Wei Wang, Prof. Wen Luo, and Prof. Shi-Lun Guo for valuable comments. We appreciate the help from Dr. Yang-Ping Shen on the simulation.